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Abstract 

The recent data from HI Collaboration on hadron multiplicity in diffractive 
DIS has been studied in the framework of perturbative QCD as a function of 
invariant diffractive mass. The formulas obtained explain the observed excess 
of particle production in diffractive DIS relative to that in DIS and e^e~ 
annihilation. It is shown that the results are sensitive to the quark-gluon 
structure of the pomeron. Namely, the data say in favour of a super-hard 
gluon distribution at the initial scale. 
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I. INTRODUCTION 



One of the most remarkable and intriguing discoveries in hadron physics in the last years 
was the observation of hard interactions (namely, partonic activity) in diffractive events by 
the UA8 Collaboration at the CERN Collider. This discovery was inspired by a seminal 
paper by Ingelman and Schlein in which such a possibility was foreseen as a consequence 
of the pomeron having an internal structure and a quark/gluon content. 

Experiments performed at HERA by ZEUS |^ and HI |^ collaborations provided further 
elements to this view through the obtainment of deep inelastic electron-proton scattering 
(DIS) events accompanied by large rapidity gaps adjacent to the proton beam direction. 
The presence of rapidity gaps in such events is interpreted as indicative that the internal 
structure of a colourless object carrying the vaccum quantum numbers, namely the pomeron, 
is being probed 0,3. 

Further experimental evidences of hard diffraction have been reported by CDF and DO 
collaborations in terms of diffractive production of W's and dijets 0,^. 

Taking together, these evidences seem to indicate that the Ingelman- Schlein (IS) pic- 
ture is right at least qualitatively. However, a serious problem arises when one checks 
this model quantitatively. For instance, the model predictions systematically overestimate 
the diffractive production rates of jets and W's [0. There are reasons to believe that the 
discrepancy between predictions and data comes from the so-called pomeron flux factor . 
In fact, a recent analysis has shown that the pomeron structure function extracted from 
HERA data by using the IS approach is strongly dependent on which expression is used for 
the flux factor. 

Recently, new data from HI Collaboration on final hadronic states in diffractive deep 



inelastic process (DDIS) of the type 

ep e'XY, (1) 

where X and Y are hadronic systems, have been presented. The system X and Y in (|I]) are 
separated by the largest rapidity gap. Y is the closest to the proton beem (My < 1.6 GeV) 



and squared momentum transfer at the pY vertex, t, is limited to |t| < 1 GeV^, while (Q^) 
is 21 ^ 27 GeV^ [|l^. Both invariant masses, Mx and My, are small compared to W, the 
centre-of-mass energy of the 7*p system. 

In particular, charged hadron multiplicity have been studied in [|T0| as a function of 
Mx in the centre-of-mass system of X. The data obtained have been compared with the 
calculations of JETSET e+e~ (which is known to reproduce well the e^e^ data) and with 
the data on hadron multiplicities in deep inelastic scattering (DIS) |Tl[]. The interesting 
observations presented in are the following: i) for Mx > 10 GeV {n)^^^^{Mx,Q'^) 
is larger than charged hadron multiplicity in DIS, {n)^'^(W,Q'^), at comparable values of 
W = (Mx)', ii) {n)^^^^{Mx,Q^) is also larger than charged hadron multiplicity in e^e~ 
annihilation, (n)^^^ (s), taken at ^/s = (Mx)- 

In the present paper we demonstrate that perturbative QCD is able to explain (qualita- 
tively and quantitavely) the more rapid growth of hadron multiplicity in DDIS. So, no mech- 
anisms besides gluon/quark jet emission with subsequent jet fragmentation into hadrons are 
needed. It is shown that the results on hadron multiplicity are very sensitive to the quark- 
gluon structure of the pomeron. An advantage of the method developed here is that it does 
not depend on the pomeron flux factor. 

The paper is organized as follows. In Section II, we present the QCD formalism for the 
description of flnal hadron states in ordinary DIS. In Section III, we extend this formalism to 
diffractive DIS and apply it to describe the HI data. Our main conclusions are summarized 
in Section IV. 



II. HADRON MULTIPLICITIES IN HARD PROCESSES IN PERTURBATIVE 

QCD 

For the time being, we cannot describe a transition of quark and gluons into hadrons in 
the framework of QCD. Nevertheless, perturbative QCD enables one to calculate an energy 
dependence of characteristics of flnal hadrons produced in hard process {e^e" annihilation 
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into hadrons, DIS, Drell-Yan process etc.). Mean hadron multiplicity, (n), is one of the 
main features of final hadronic states. Hadron multiplicity in e^e~ annihilation has been 
studied in a number of papers. The result looks like (see, for instance, ||12||) 



{nY 



a(lns)^exp {cy/his^ 



(2) 



where ^/s is a total c.m.s. energy of colliding leptons. As one can see, {nY ^ {s) rises more 
rapidly than In s, although more slowly than any power of s. Expression describes well 
the available data. 

Hadron multiplicity in DIS was calculated first in the framework of perturbative QCD 
in Refs. (see also 10). It was shown that average multiplicity in lepton scattering off 
parton 



eq{g) e'X, 



(3) 



(n 



pis 

'q/g 



(W^,g2), is related to {nY^ (s), the average multiplicity in e^e annihilation, taken 



at y/s = W (up to small NLO corrections which descrease in W and Q"^) [|T3|. 

In a case when quark distribution dominates (say, at x ~ 1), the relation between {n)^^^ 
and (nY^'^ looks like (up to small NLO corrections which decrease in W and Q"^) ITi 



DIS 



(4) 



If we consider small x, we have to account for the gluon distribution and the result is of 



the following form ||15 



DIS 



{n)^ iW,Q') = {nY^^ (W) 



1 H — —s 1 e 

2Cf \ 2 



(5) 



where Ca = Nc, Cp = {N^ — i)/'2Nc and Nc is a number of colours. 
The quantity e is defined via gluon distribution (see for details) 



lnDf(e,x) 



(6) 



V 2-nCA di 

^ being the QCD-evolution parameter ^(W^^) = [dk"^ / k'^){as{k'^) /27t) . Let us notice that 
e (^ does not depend on a type of a target, and e is completely defined by the evolution of 
D9 in ^. 



Starting from the well-known expression for at small x [16 

1 



ln(l/x 



-f Ii(2f ) exp(— (i^), 



(7) 



where v = y 4A^cC d = (1/6) (llA^c + "^Nf/N^), Nf is a number of flavours, we get 

e = e{W^,x) -- 



ln(l/x) 



d 



(8) 



At ln(l/x) ^ 1 (that is omitting the second term in RHS of (§)) we come to the asymptotical 



expression for e from Ref. |15 



Using (^) and (|), we get in a general case 



DIS . 



{n) {W,Q') = {nr^-iW) 



A + ( 1 + ^e{W^ e{W\ x) ) ) (1 - A) 



(9) 



where A/(l — A) defines the quark/gluon ratio inside the nucleon. As was shown in e 
(§) is a leading correction to (n) which rises with the decreasing of x. 

All mentioned above is related to the multiple production of the hadrons in DIS of lepton 
off the parton (|]). It is a subprocess of the process of lepton deep inelastic scattering off the 
nucleon 



ep 



e'X. 



(10) 



According to Refs. [|^,0,|1^, the corresponding formula for mean hadron multiplicity in 
DIS looks like 



^ D/5 

{nf{W, Q') = (n) {W^ff, Q') + n,{M',). 

DIS 



(11) 



The quantitiy (n) has been defined above (see (|^)). It depends on the effective energy, 
Weff, available for hadron production in the parton subprocess (|^): 

kW\ (12) 

The average multiplicity of nucleon remnant fragments, ng, is a function of its invariant 
mass Mq = (1 — z){Ql/z + M^), where M is a nucleon mass. 

The efficiency factor k in (|12|) was estimated in to be (k) ~ 0.15-^0.20 at present en- 
ergies. That is why {n)^^^(W, Q'^) is less than {nY'^''~ {W) and the growth of {n)^^^{W, Q"^) 
in W at fixed is delayed in comparison with that of (n)*^^^ (W^)- 



III. AVERAGE HADRON MULTIPLICITY IN DIFFRACTIVE DIS 



Hadronic system X in diffractive DIS (0) is produced as a result of the virtual photon- 
pomeron interaction: 

eP ^ e'X. (13) 
The kinematical variables usually used to describe DDIS (in addition to W and Q"^) are 

We start from the fact that pomeron has quark-gluon structure. It means that hadron 
production in process (|T3) is similar to that in parton subprocess of DIS (Bf). 



Recently a factorization theorem for DDIS has been proven |T^. Structure functions of 
DDIS coincide with structure functions of DIS and have the same dependence on a factor- 
ization scale. As a result, distribution functions of quark and gluons inside the pomeron, 
-Dp^\ must obey standard DGLAP evolution equations at high |]20[| . 

So, we can conclude that a hadron multiplicity of system X in DDIS, {n)^^^^ , is given 
by expression (Rl), in which W is replaced by Mx, while x is replaced by /?. Namely, we get 



(n)^^^^(Mx,Q') = (n)'^^^-(Mx)x 

(15) 



X 



Ap + ( 1 + ^e{Ml e{Ml, /3))) (1 - Ap) 



2Cf 

Here e(M|^,/3) is given by formula (§) and Ap/(1 — Ap) defines quark/gluon ratio inside 
the pomeron. As both Dp and -Dp obey DGLAP equations, the quantity e(M|^,/3) and, 
consequently, hadronic multiplicity in DDIS are sensitive to the quark-gluon structure of 
the pomeron at starting scale Qq. 
As folUows from (|ll|) and (|15|), 

(n)^^^^(Mx, Q') > {n)^''{Mx, Q') (16) 

because {n)^^^^{Mx, Q"^) is defined by (n)^"^"' (Mx), while (n)^^'^(Mx, Q"^) is expressed via 
l^YiY*^ [M'^^), the quantity M'^^ being much smaller than Mx by factor k (see (|12D). 



For Mx < 29 GeV we have M*^^ < 5 ^ 6 GeV. It is known that in the region W <b^Q 
GeV function (n)^^'^ {W) rises logarithmically (~ Inl^) that is slower than (|^). This results 
in more rapid growth of {n)^^^^{Mx,Q'^) in Mx in comparison with {n)^^^{Mx,Q'^)- 

We conclude from formula (|15D that {n)^^^'^{Mx, Q^) should exceed {nY^^^' {Mx)- More- 
over, the ratio {n)^^-^^{Mx,Q'^)/ {n)'^^^ (Mx) has to grow in Mx at fixed (that corre- 
sponds to the rise in at fixed Q^). All said above is in good qualitative agreement with 
the data from HI Collaboration |T0[. 

In Fig. 1 we show the result of the fits of the HI data by using our formula (^) (solid 
curve). In order to obtain this result, we proceeded as follows. For the quantity Ap which 
enters expression (pTsD we used 

where and are respectively the gluon and the singlet quark distributions inside 
the pomeron that, as mentioned above, obey DGLAP evolution equations For the 

distributions at the initial scale Ql = 4 GeV^ the following forms were employed: 

DUP,Ql)=ai P{l-P), 

DUf3,Ql)=h(3''%l-f3f\ (18) 

Thus, for the quark distribution we fixed an initial hard profile leaving free the normalization 
parameter, while for the the gluon distribution we have left all parameters free without 
imposing any sum rule. In order to perform the fit, other elements are needed. In Eq. (|15|), 
for {nY"^^ {Mx) we have used the parametrization 

(n)^^^"(Mx) = 2.392 + 0.024 ln(-J^) + 0.913 ln2(-j:) (19) 



with Mo = 1 GeV ^ taken from Besides that, the experimental data shown in Fig.l are 
given in terms of average values of Mx, and that do not obey strictly the kinematical 
relation (0). In order to be faithful to data, we employed in our fit the parametrization 



with a = 1.63, b = 0.165 GeV-\ and c = 2.202. 

The dotted (dashed) hne in Fig. 1 corresponds to a pure quark (gluon) content of the 
pomeron. The sohd hne, which gives the fit of the HI data, is obtained 

with parameters ai = 2.400, h = 3.600, 62 = 5.279, and 63 = 0.204 in Eqs. (|T8|), which 
are evolved to the respective Q^- value corresponding to each experimental point. 

The initial distributions and their evolution with are shown in Fig. 2. As can be seen, 
a super-hard profile (peaked at /3 ~ 1) was obtained for the gluon distribution at the initial 
scale, in qualitative agreeement with HI analysis p2| . 

The Q^-evolution of the normalized fractions of the pomeron momentum carried by 
quarks and gluons is presented in Fig. 3. It is shown that quarks are slightly predominant 
over gluons at the initial scale, but that this relation reverses as evolves. In this case, our 
results do not follow those obtained by HI Collaboration since their analysis [|10[ favour a 
fit in which predominate gluons carrying > 80% of pomeron momentum at the initial scale. 
It must be said, however, that the results of Fig. 3 are consistent with limits established in 
other experiments (see, for instance, [Q). 



IV. CONCLUSION 

We have presented in this paper a description of the hadron multiplicity obtained in 
diffractive DIS and recently reported by the HI Collaboration ||10[. This description was 
derived from a formalism previously developed for ordinary DIS within the framework of 
the perturbative QCD. 

The formula obtained enables us to explain the observed excess of hadron multiplicity in 
diffractive DIS, {n)^^^^{Mx,Q'^), relative to those in DIS and e~^e~ annihilation taken at 
W = (Mx) and ^/s = {Mx), correspondingly. The more rapid growth of {n)^^^^{Mx,Q'^) 
is also understood. 

It was shown (Fig.l) that neither a pure quark nor a pure gluon content of the pomeron 
satisfy the data, but that a mixing of both components in approximately equal shares is able 
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to provide a good description of the experimental results. The pomeron structure function 
that comes out from the present analysis consists of a hard distribution for the quark singlet 
and a super-hard distribution for gluons at the initial scale of evolution, in agreement with 



what has been reported lately in the literature [Q. 

This result is remarkably good if one considers that it was obtained from a very small 
data set (only 7 data points), covering a limited (3 range (0.03 < /3 < 0.43) for an equally 
restricted Q^.^ange (22 < < 27 GeV^). 

Concluding, the most important result presented here is the theoretical framework sum- 
marized by Eq.(|l^), that has made it possible to explain anomalous HI data on hadron 
multiplicity in DDIS. It also enabled us to extract information about the pomeron structure 
function from such a limited data set, without the usual complications and ambiguities with 
flux factor normalization and xp-dependence. 

This conclusion points out to the need of more and more precise DDIS multiplicity data, 
taken at extended kinematical ranges. Such a possibility would improve a lot the analytical 
capacity of the scheme presented here. 
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FIGURE CAPTIONS 



Figure 1: Description of the average multiplicity of charged hadrons produced in diffrac- 
tive DIS experiments. Data obtained by the HI Collaboration ||10[. The theoretical results 
are obtained from Eq.([T5|) with different assumptions for the quark/gluon content of the 
pomeron: pure quark (dotted curve), pure gluon (dashed curve) and both components (solid 
curve) . 

Figure 2: Q"^ evolution of the quark singlet (a) and gluon (b) distributions obtained by the 
fitting procedure described in the text. 

Figure 3: evolution of the fractions of the pomeron momentum carried by quarks and 
gluons as predicted from the parametrization resulted from the fit to the HI data. 



12 



Figure 1 
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